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ABSTRACT 


Finite  amplitude  acoustic  standing  waves  in  a  rectangular 
air-filled  cavity  with  various  wedge-shape  boundary  perturbations 
were  studied  both  experimentally  and  theoretically.  The  experi¬ 
mental  results  show  that  geometrical  perturbations  alter  the 
finite-amplitude  behavior  of  the  cavity  and  that  the  nature  of 
these  changes  are  in  qualitative  agreement  with  the  predictions  of 
the  theory.  However,  quantitative  agreement  was  not  observed, 
possibly  because  the  perturbation  chosen  did  not  satisfy  all  the 
assumptions  of  theory. 
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I.  INTRODUCTION 


Coppens  and  Sanders  Cl,  2]  developed  a  non-linear  acoustic 
model  with  the  dissipative  term  phenomenologically  describing  the 
viscous  and  thermal  energy  losses  actually  encountered  at  the 
walls  of  a  rectangular  rigid  cavity.  Several  researchers  [3,4] 
have  examined  the  problems  and  experimental  results  are  in 
excellent  agreement  with  the  theory  except  when  a  degeneracy 
exists  in  the  cavity. 

The  purpose  of  this  research  was  to  examine  the 
finite-amplitude  behavior  within  a  cavity  for  which  perturbation 
effects  can  be  accurately  measured  and  to  compare  the  results  with 
the  predictions  of  the  model. 


II.  BACKGROUND  AND  THEORY 


In  1975  Copperas  and  Sanders  C 2]  formulated  a  perturbation 
expansion  for  the  non-linear  acoustic  wave  equation  with  a  dissi¬ 
pative  term  describing  the  measured  absorption  properties  and  the 
measured  resonance  frequencies  for  standing  waves  within  a  real, 
fluid  filled,  rigid-walled  cavity.  This  model  predicts  that,  when 
a  near-degeneracy  exists,  geometrical  perturbations  provide  a 
mechanism  whereby  a  nearly  degenerate  mode  can  affect  the 
finite-amplitude  behavior  of  the  cavity. 

The  non-linear  wave  equation  of  the  viscous  fluid  [2]  is, 

(co2Q2+3</3t)  p/pQco2«32/3t2[(u/co)2+»s(Y-l)  (P/P0c02)2] 

where  cQ2  m  (3p/3p)  (adiabatic},  pQ  is  the  equilibrium 
density  of  the  fluid,  p  the  acoustic  pressure,  u  the  magnitude  of 
the  acoustic  velocity,  y  the  ratio  of  heat  capacities, and <£  an 
operator  describing  the  physical  processes  for  absorption  and 
dispersion.  The  term  on  the  right  can  be  interpreted  as  a 
distribution  of  virtual  sources  created  by  the  self-interaction 
of  the  standing  waves. 

Pressure  standing  waves  in  a  rectangular  cavity  of  dimension 
L*,  Ly,  Lz  have  the  form 

p  cos  kxx  cos  kyy  cos  kzz  cos  wt 

where 

kx  »  jx/Lx  j  “  0,1,2,... 

Ky  *  (ir/Ly  I  *  0,1,2,... 

Kx  *  raiT /bz  a  *  0,1,2,... 

and  j+f+mjK) 


If  the  cavity  ia  being  driven  near  resonance,  the  contributions  of 
non-resonant  terms  are  negligible  with  respect  to  the  resonant 
terms;  the  acoustic  field  has  the  form 


p*  £pn 

itfO 


pr/poco^®n  cosniVc  oosnkzz  sin(nwt+$n) 


(2) 


M  is  the  Mach  number  |  u/cj  ,  R]_»l,  and  Rn  is  the  relative 
amplitude  of  the  n-th  standing  wave.  By  substituting  (2)  into 
(1),  a  set  of  coupled,  non-linear,  transcendental  equations  is 
obtained , 

\  fS?  <w  vj 

(3) 

-  ^  fflVi1] 

where  H  »  1  for  axial,  2  for  tangential,  and  3  for  oblique 
standing  waves. 

If  a  geometrically-perfect  rectangular  cavity  is  driven  at 
frequencies  ne ar  the  resonance  frequency  of  the  (0,1,0)  standing 
wave,  only  the  family  members  (0,n,0)  will  contribute 
significantly  to  the  finite  amplitude  behavior. 

In  this  research  we  were  interested  in  this  behavior  when  a 
degeneracy  exists  between  the  (0,2,0)  and  (1,0,0)  modes  and  the 
cavity  is  not  a  perfect  right  parallelepiped. 

If  we  define  a  perturbation  parameter  £  as  a  dimensionless 


measure  of  the  magnitude  of  any  irregularities  on  the  cavity 


surfaces  compared  to  the  effective  dimensions  of  the  cavity,  the 
effect  of  the  perturbation  on  the  standing  wave  Pn  that  would 
exist  in  the  ideal  cavity  can  be  expressed  as  a  small  correction 
p1  so  that  the  true  pressure  field  pn'  is 
Pn‘  "  Pn  +  £P' 

For  the  case  of  interest,  n  denotes  the  (0,2,0)  standing  wave  and 
p*  is  the  perturbation-generated  (1,0,0)  standing  wave  of  nearly 
identical  resonance  frequency. 


III.  SAMPLE  CALCULATION  OF  THE  PRESSURE  DISTRIBUTION  WITH  A 
WEDGE  PERTURBATION 

The  pressure  distribution  of  the  (0,2,0)  wave  has  the  form 

P2  ■  p  COS  (2iry/Ly)  cos(2a>t+92)  (1) 

where  P  and  92  are  constant  determined  by  the  driving 
conditions . 

For  this  example,  let  the  equation  for  the  perturbed  boundary 
(Fig.  o)  be 

X  “  LX  (1“(A/Lx)  C (4/V  (yw3V4)  ]  KHy-3Ly/4)-0(y-I^)  ]} 

-  Lx  [H€f(y,z)]  {2) 

where 


£H4/lx 

f(y,2)  =“  -(4Ay)  (y-3Ly/4)  [0(7-31^4) -Uty-I^U 


U(y) 


0  if  y  <  0 
1  if  y  >  0 


Differentiating  p2  and  f(y,z)  with  respect  to  y  gives 


3P2/3Y  *  -(2irP/Ly)  sin(27ry/Ly)  cos(2cot+92)  __  ' 

and 

3f/3y  =  -(4/Ly) [U(y-3Ly/4)-U(y-Ly) ] 

-(4/L  ) (y-3Ly/4) [5 (y-3Ly/4) -5 (y-Ly) ]  (3) 

and  5 (y)  is  the  Dirac  delta  function. 

Now,  it  can  be  shown  by  standard  perturbation  analysis  [4]  that 
the  perturbation  correction  p'  must  satisfy  the  boundary 
condition 

(3pV3x)L  -  Lx(3f/3y)  (3p2/3y)L 

X  X 

Substituting,  we  obtain 

(3p'/3x)  *  A  sin(2iry/L)  cos  (2cot+9  ,)  (u  (y-3L  /4)  -U  Cy-L  ) 

^x  Y  *  y  y 

+  (y-3Ly/4)  [5 (y-3L  /4)-« (y-L  ) ]>  (4) 

where 

A  -  8irPLx/Ly2 

Equation  (4)  can  be  expanded  as  a  Fourier  series, 

00 

(3p'/3x)T  *  A  cos(2a)t+9,)  Z  [A_  cos  (jmry/L  )  +b_  sin  {my/1‘,  A  J  C5) 

"v  i  m=n  m  y  m  y 
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Only  Aq  is  nesdsd  to  determine  ths  first  ordsr  corrsction  term 
dus  to  ths  degeneracy  betwesn  ths  (0,2,0)  and  (1,0,0)  waves: 


AQ  -  1/2* 

Ths  psrturbation  corrsction  p*  is 

p'  *  (Ac02/irLx4o>2)  Qi00sinT100  coa  t^/^)  cos  (2a>t+92+T100)  (6) 

whsrs  ths  phass  angle  T  is  given  by 

tanTioo  "  W^W2^*1 

The  pressure  at  the  microphone  position  x*Lx  can  be  written  as 
follows  (with  2«/c0»4r/Ly) , 

pmic  "  p2  +  ep* 

*  P(cos(2tot+62)  +  (c/2tt  2)  Q1O(j®inTiooCO8(2“t+02+TlOOi]  (7) 

The  total  pressure  amplitude  at  ths  microphone  will  be 

I  Pmid  *  P01+B  cos  Tioo)2  ”  (B  ain  T100)2  ^81 

where 

B  -  (e/2x2)  qiqo  8in  rl00  C9a) 

B  can  also  be  expressed  in  terms  of  Ag  for  a  wedge  anywhere  on 
the  wall  x*Lx  by 

B  -  (e  Aq/x)  qi00  sin  ti00  (9b) 
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where  \  is  calculated  for  any  arbitrary  position  and  wedge 
dimensions  by  tbs  same  method  as  developed  above. 

Cln  reference  [$],  the  author  made  a  error  in  formulating  Equation 
(2);  he  used  y  instead  of  (y-aL*)  and  got  an  incorrect  value  of 
\j*  3 
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DEFINITION  OF  SOME  PARAMETERS 


Several  parameters  will  be  defined  for  the  purpose  of 
simplifying  the  mathematical  formulation  and  elucidating  the 
physical  content  of  the  equations. 

A.  STRENGTH  PARAMETERS:  S 

S  -  MgQ^ 

where  ?  «  (1/2)  (1+y)  and  is  the  quality  factor  of  the  driven 
(0,1,0)  standing  wave.  The  strength  parameter  characterizes  the 
strength  of  the  finite-amplitude  interaction.  It  is  interesting 
to  note  that  S  is  one  half  the  Goldberg  number.  Since  the 
microphone  sensitivity  Sm  (obtained  with  a  B&K  4220  pistonphone) 
is  known, 

S  ■  MgQi  »  7.07  x  10“3  VjQ! 

where 

m  «  yfl  v1/pc02sm,  8  *  1.2  for  air, 

P  ■  1.293  kg/m3,  cQ  ■  345  m/s 

Vx  «  RMS  output  voltage  of  the  first  harmonic  component, 

and 

Ql  is  the  quality  factor  of  the  driven  (0,1,0)  standing 
wave . 

B.  FREQUENCY  PARAMETER:  Fn 

F„(f)  -  QnCl-(f„/f)2] 

where  f  is  the  driving  frequency,  and  fn  and  Qn  are  the 
resonance  frequency  and  quality  factor  for  the  (0,n,0)  wave. 

If  Fn<l,  the  nth  harmonic  of  the  driving  frequency  lies 


within  the  half-power  frequencies  of  the  resonance  curve  for  the 
(0,n,0)  wave. 


For  reasonably  large  values  of  Q^, 

Fi  (f)  •  2  QX 

C.  HARMON 1C ITY  COEFFICIENT;  E(n) 

E(n)  ■  {fn-nfi)/nfi 

E(n)  characterizes  how  well  the  modes  of  a  given  family  are 
tuned  (harmonic).  If  |  E(n)|  <4,  the  corresponding  harmonic  will  be 
strongly  excited. 


LO 


V.  APPARATUS 


A.  THE  RECTANGULAR  CAVITY 

The  cavity  (Fig.l),  constructed  from  0.75-in.  aluminum,  has 
interior  dimensions  12.00  in.  long,  2.50  in.  high,  and  a  width 
that  can  be  varied  between  5.50  in.  to  7.00  in.  in  0.25  in. 
increments.  All  joints  were  sealed  with  a  thin  layer  of  silicon 

grease . 

Figure  2  shows  the  pressure  distribution  for  several  of  the 
lower  modes  of  this  cavity,  and  Table  0  presents  the  theoretical 
eigen  frequencies  calculated  for  a  cavity  12  x  6  x  2.5  in.  Note 
that  the  (0,2,0)  and  (1,0,0,)  modes  are  predicted  to  be 
degenerate.  Experimentally,  the  resonance  frequency  of  the 
(0,2,0)  standing  wave  was  about  3  Hz  higher  than  that  of  the 
(1,0,0)  standing  wave.  The  various  configurations  of  wedge- 
perturbations  are  shown  in  Fig.  3. 

The  source  piston  is  set  flush  with  the  bottom  of  the  cavity 
as  near  to  the  wall  y*0  as  practical,  and  halfway  between  the  wall 
at  x«0  and  x*Lx.  In  this  position  it  can  efficiently  excite  the 
#(0,n,0)  family  of  waves  without  appreciably  exciting  the  (1,0,0) 
wave.  To  determine  the  proportion  of  the  (1,0,0)  mode,  an 
auxilliary  driver,  an  (ID-30),  can  be  inserted  at  Position  A. 

A  microphone  at  Position  B  will  sense  the  pressure  of  all 
standing  waves.  If  the  microphone  is  placed  at  Position  A,  it 
senses  the  (1,0,0)  wave  with  minimum  contamination  from  the 
(0,2,0)  wave  and  at  Position  C  it  senses  the  (0,2,0)  wave  with 
minimum  contamination  from  the  (1,0,0)  wave. 
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B.  APPARATUS 


A  block  diagram  of  the  apparatus  Is  shown  in  Fig. 4.  A  GR 
1161 -A  coherent  decade  frequency  synthesizer  is  used  to  produce  a 
driving  signal  precise  to  within  ±0.001  Hz.  This  signal  is 
applied  to  a  2120  MB  power  amplifier  which,  in  turn,  drives  the 
shaker.  The  motion  of  the  piston  was  continuously  sensed  by  an 
Cndevco  Model  2215  Accelerometer  mounted  within  the  piston,  and 
the  output  of  the  accelerometer  was  observed  on  a  Model  130BR  HP 
oscilloscope  and  measured  on  a  HP  400D  Voltmeter.  A  Schlumberger 
spectrum  analyzer  was  used  to  measure  the  harmonic  distortion  in 
the  piston  motion. 

The  output  of  the  B4K  1/4-in.  microphone  (with  matching 
preamplifier  B&K  2810)  was  fed  into  three  devices;  (1)  an  HP  400D 
VTVM  to  measure  overall  voltage  level,  (2)  a  Schlumberger  spectrum 
analyzer  to  display  the  spectrum  of  the  waveform,  and  (3)  two  HP 
302A  wave  analyzers  to  measure  the  amplitudes  of  the  first  two 
harmonics  of  the  pressure  waveform. 
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VI.  DATA  COLLECTION  PROCEDURE 


Since  the  resonance  frequencies  of  the  cavity  were  observed 
to  vary  with  time,  the  system  was  allowed  to  warm-up  for  at  least 
one  hour  prior  to  data  collection.  The  piston  was  then  driven  at 
the  maximum  amplitude  to  be  expected  during  the  run,  and  the 
harmonic  content  of  the  accelerometer  output  was  analyzed  and  the 
piston  was  adjusted  until  the  second  harmonic  of  the  accelerometer 
output  was  at  least  50  dB  below  that  of  the  fundamental.  Figure  5 
shows  typical  results  for  the  percent  second  harmonic  in  the 
acceleration  with  the  piston  driven  at  a  rather  large  amplitude. 
Finite-anplitude  measurements  were  limited  to  those  frequencies 
for  which  V2/Vi<0.01. 

There  are  three  steps  necessary  for  collecting  accurate  data; 
(1)  pre-run  infinitesimal-amplitude  measurements,  (2)  finite- 
amplitude  measurements,  and  (3)  post-run  infinitesimal-amplitude 
measurements . 

A.  PRE-RUN  AND  POST-RUN  INFINITESIMAL-AMPLITUDE  MEASUREMENTS 

During  these  measurements,  the  piston  was  driven  with  the 
accelerometer  output  less  than  0.1  V.  Observation  of  spectra  of 
the  pressure  waveforms  obtained  at  these  low  amplitudes  showed 
that  the  amplitudes  of  all  overtones  were  at  least  60  dB  below 
that  of  the  fundamental.  The  resonant  frequencies  fn  and  the 
quality  factors  Qn  for  the  first  few  (0,n,0)  waves  were 
determined  from 
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v 

y:  ■ 


fn  -  <fu+f*)/2 


On  -  *n/<*u-**> 

where  fu  and  f^  ara  tha  uppar  and  lower  ha If -power  frequencies. 

The  resonance  frequency  and  Q  of  the  (1,0,0)  wave  were  determined 
by  the  same  procedure  but  with  the  cavity  driven  at  port  A.. 

The  harmonicity  coefficients  were  then  calculated  from 

En  *  ( f n“»f i ) /nf i ,  for  (0,n,0)  wave 
•  (fl00~f2)/f2*  for  the  (1,0,0)  wave 

where  fioo  i®  the  resonance  frequency  of  the  (1,0,0)  wave. 

The  time  at  which  each  fu  and  were  taken  was  recorded. 

B.  THE  FINITE  AMPLITUDE  MEASUREMENTS 

Throughout  this  portion  of  the  run,  the  strength  parameter  S 
was  maintained  constant  by  adjusting  the  driving  voltage  applied 
to  the  piston;  with  the  microphone  in  Position  B,  the  HP  302A  wave 
analyzer  was  set  to  the  driving  frequency  and  the  driving  voltage 
adjusted  to  keep  the  amplitude  of  the  fundamental  of  the  pressure 
waveform  constant  as  the  frequency  was  changed.  The 
frequency  was  increased  0.5  Hz  steps  through  ±6  Hz  about  the 
resonance  frequency  of  the  fundamental.  At  each  driving 
frequency,  the  harmonic  content  of  the  microphone  output  (usually 
up  to  the  fourth  harmonic)  was  measured  with  the  second  HP  302A 
wave  analyzer  set  on  AFC  mode.  The  time  of  each  measurement  was 
recorded . 
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To  determine  the  frequency  parameter  at  tha  tiro*  the  harmonic 
content  was  measured,  fn  at  the  instant  of  the  measurement  of 
Vn  was  estimated  by  interpolation  of  the  fn  found  in  the  pre- 
and  post-run  procedures.  Figure  6  shows  the  drift  in  fn  is 
approximately  linear  in  time.  The  results  were  presented  as 
Vn/Vi  vs  Pj_(f)  for  a  given  strength  parameter  and 
perturbation. 


VII.  RESULTS 

The  results  aura  presantad  in  graphical  form  in  Appandix  A  and 
in  tabular  form  in  Appandix  B. 

Figure  7  shows  the  excellent  agreement  between  theory  and 
experiment  obtained  when  there  are  no  perturbations  so  that  only 
one  family  of  waves  is  excited.  Some  discrepancies  were  always 
observed  when  the  driving  frequency  was  far  from  the  fundamental 
resonance  frequency  (for  |  FjJ  >2).  In  this  region,  the  piston 
had  to  be  driven  hard  to  keep  Vj_  constant,  causing  significant 
second  harmonic  to  appear  in  the  piston  waveform  (as  shown  by 
Fig.  5 ) ,  thereby  introducing  a  linearly-generated  second  harmonic 
into  the  cavity. 

For  all  cases  studied,  the  pressure  was  calculated  following 
the  methodology  of  Sect.  III. 

For  the  boundary  perturbation  of  Fig.  3,  the  results  are 
plotted  in  Figs.  8  through  15.  The  thin  line  is  the  theoretical 
prediction  in  the  absence  of  any  perturbation  and  the  thick  line 
includes  the  perturbation  correction.  The  circles  are  the 
experimental  results.  Along  the  edge  of  each  figure  is  the 
information  about  the  cavity  configuration. 

A.  THE  WEDGE  AT  THE  CORNER  OF  THE  CAVITY 

Figures  8  through  11  show  the  results  for  a  wedge  at  the 
corner  of  the  cavity.  The  wedge  used  are  those  of  Fig.  3(a)-(d) 
respectively. 

B.  THE  WEDGE  AT  THE  CENTER  OF  THE  LONG  WALL 

Figures  12  through  15  show  the  results  for  the  same  wedges 
[Fig.  3(e) -(h)]  but  now  located  at  the  center  of  the  long  wall. 


For  Fig.  12  and  13,  AQ  -  0.177  with  &  for  Fig.  13  twice  that  for 
Fig.  12.  For  Fig.  14  and  15,  Aq  -  0.500  with  A  the  same  as  for 


VIII.  CONCLUSIONS 


A.  The  experimental  apparatus  and  procedures  are  capable  of 
providing  data  sufficiently  precise  to  verify  the  prediction  of 
the  theory  in  the  absence  of  geometrical  perturbation  if  the  range 
of  driving  frequencies  is  restricted  so  that  the  absolute  value  of 
the  frequency  parameter  is  less  than  2. 

B.  To  verify  the  correction  for  a  geometrical  perturbation, 
it  is  useful  (1)  to  have  a  sufficiently  large  correction  to  the 
unperturbed  prediction  and  (2)  that  the  sign  of  the  correction 
change  for  an  absolute  value  of  the  frequency  parameter  less  than 
2. 

C.  The  experimental  results  show  that  geometrical 
perturbation  alters  the  finite  amplitude  behavior  of  the  cavity, 
and  that  the  nature  these  changes  are  in  qualitative  agreement 
with  the  predictions  of  the  theory.  However,  quantitative 
agreement  was  Jt  observed. 

D.  Sources  of  the  difficulties  in  obtaining  good  agreement 
might  be  (1)  the  inability  to  experimentally  satisfy  conditions  of 
B  above,  and  (2)  the  higher  order  terms  neglected  in  developing 
the  theory  may  not  all  be  small.  For  example,  terms  of  order 
higher  than  first  order  in  Eq.  (  ♦)  of  Sect.  II  and  terms  of  order 
higher  than  Aq  and  all  b^  in  Eq.  (5)  of  Sect.  IV. 


Fig.  1.  Adjustable  Rectangular  Cavity 


Block  Diagram 
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Fig.  13 


Fig.  15 


TABLE  0 


Mode 

Frequency  (Hz) 

Mode 

Frequency  (Hz) 
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565 
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020 

1131 

140 

2529.0 

030 
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3032.0 

100 
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2757.5 
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Oil 

2926.4 

TABLE  I  (Fig.  7) 
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time 

fu 
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f 

r 

Q 

E 

FP 
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020 
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100 

Vl* 

0 

3 

6 

9 

lldB 

565.77 

1131.84 

1698.27 

1131.87 

time  f 
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-2.12 
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20 

562.5 
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0.0380 

0.0025 

-1.43 

21 
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-85.0 
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0.0052 
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-0.05 

25 
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0.0200 

0.0067 

0.29 

26 

565.5 
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0.0081 

0.64 

27 
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0.0610 
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31 
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32 
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-44.1 
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33 
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3.05 
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0.0013 
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-67.2 
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41 
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TABLE  II  (Fig.  8) 
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189.81 

0 

020 

3 

1132.37 

1127.62 

1129.99 

238.25 

-0. 31xl0-2 

030 

6 

1701.89 

1696.79 

1699.34 

333.14 

-0.57x10-3 

100 

9 

1132.44 

1127.71 

1130.07 

238.76 

0.69xl0~4 

VjS-lldB  time  f 

V2 

V  V 

3  4 

V*! 

V*1  P 

12 

561.3 

-44.3 

-69.9 

-71.3 

13 

561.8 

-43.2 

-69.2 

-72.2 

14 

562.3 

-41.7 

-68.3 

-74.0 

15 

562.8 

-40.2 

-66.8 

-80.0 

16 

563.3 

-38.5 

-65.2 

17 

563.8 

-36.4 

-62.9 

18 

564.3 

-34.3 

-60.3 

19 

564.8 

-32.8 

-57.2 

20 

565.3 

-32.5 

-55.0 

-72.0 

21 

565.8 

-33.5 

-53.5 

-67.8 

22 

566.3 

-35.1 

-51.8 

-62.5 

23 

566.8 

-36.9 

-50.9 

-60.9 

24 

567.3 

-37.9 

-51.2 

-60.8 

25 

567.8 

-38.6 

-52.8 

-64.6 

26 

568.3 

-39.2 

-55.3 

-69.2 

27 

568.8 

-40.3 

-58.7 

-75.0 

28 

569.3 

-41.5 

-61.7 

-78.5 

29 

569.8 

-42.4 

-64.3 

— — 

30 

570.3 

-44.3 
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0.0024 
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0.0032 

0.15 
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0.49 

0.0072 

0.83 
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1.86 
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TABLE  IV  (Fig.  10) 
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563.20 

564.70 

189.30 

0 
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3 
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1124.95 
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287.70 

-0.2x10 
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6 
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-0.18xl0~3 
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9 
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time  f 

V2 
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V*1 

V*!  FP 

12 

558.7 

-42.8 

-69.8 

0.0257 

0.0011 

13 

559.2 

-41.8 

-68.7 

0.0288 

0.0013 

14 

559.7 

-40.8 

-67.9 

0.0324 

0.0014 

15 

560.2 

-39.7 

-66.9 

0.0367 

0.0016 

16 

560.7 

-38.3 

-65.5 

0.0434 

0.0019 

17 

561.2 

-36.7 

-63.5 

0.0519 

0. 0024 

18 

561.7 

-34.9 

-61.0 

0.0638 

0.0032 

19 

562.2 

-33.0 

-57.8 

-76.3 

0.9794 

0.0046 

20 

562.7 

-30.7 

-53.7 

-70.7 

0.1041 

0.0074 

0.0010 

21 

563.2 

-28.9 

-50.2 

-65.3 

0.1274 

0.0110 

0.0019 

22 

563.7 

-28.5 

-47.5 

-61.4 

0.1334 

0.0150 

0.0030 

23 

564.2 

-30.4 

047.2 

-59.0 

0.1078 

0.0155 

0.0040 

24 

564.7 

-32.2 

-48.5 

-58.5 

0.0873 

0.0133 

0.0042 

25 

565.2 

-33.7 

-52.2 

-62.7 

0.0733 

0.0087 

0.0026 

26 

565.7 

-35.6 

-56.4 

-70.07 

0.0592 

0.0054 

0.0010 

27 

566.2 

-37.0 

-58.4 

-74.7 

0.0501 

0.0043 

28 

566.7 

-39.6 

-59.7 

-77.2 

0.0314 

0.0037 

29 

567.2 

-42.1 

-61.5 

0.0279 

0.0030 

30 

567.7 

-43.7 

-63.7 

0.0232 

0.0023 

31 

568.2 

-44.5 

-65.5 

0.0221 

0.0019 

32 

568.7 

-45.1 

-66.9 

0.0197 

0.0016 

33 

569.2 

-45.6 

-68.0 

— 

0.0186 

0.0014 

34 

569.7 

-46.1 

-69.0 

— 

0.0176 

0.0013 

35 
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-46.6 

-69.8 

- — ■ 

0.0167 

0.0012 

36 
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-47.0 

-70.4 

0.0158 

0.0011 

37 
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-47.2 

-70.6 

— 

0.0155 

0.0011 

38 
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-47.7 

-71.3 

0.0146 

0.0010 
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-48.6 

-63.9 

— 
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0.0023 
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-3.71 

-3.37 

-3.04 

-2.71 

-2.38 

-2.04 

-1.71 

-1.38 

-1.04 

-0.71 

-0.38 

-0.05 

0.29 

0.62 

0.95 

1.29 

1.62 

1.95 

2.29 

2.62 

2.95 

3.29 

3.62 

3.95 
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563.86 
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3 
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1123.28 

1125.23 

288.08 

-0.49xl0“2  i 
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9 
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i 
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562.9 
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564.9 
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-49.9 

-50.0 

-50.3 

-50.4 

-50.4 

-51.1 


-66.9 

-65.3 

-63.9 

-61.8 

-59.1 

-55.7 

-55.6 

-50.8 

-50.0 

-49.7 

-50.9 

-54.3 

-59.3 

-65.0 

-69.7 

-72.2 

-74.5 

-75.4 

-73.5 

-71.5 

-70.9 

-71.1 

-71.6 

-71.8 

-72.3 

-72.4 

-72.4 

-71.6 


-78.0 

-72.7 

-68.2 

-65.3 

-62.5 

-60.3 

-61.4 

-67.0 

-75.2 


566.95  563.96  565.46 
1127.31  1123.40  1125.37 
1695.54  1690.73  1093.14 
1126.98  1123.12  1125.05 


0.0376 

0.045 

0.054 

0.0672 

0.0856 

0.109 

0.131 

0.127 

0.101 

0.081 

0.065 

0.055 

0.048 

0.043 

0.039 

0.037 

0.033 

0.030 

0.030 

0.022 

0.015 

0.012 

0.011 

0.011 

0.011 

0.011 

0.011 

0.010 


189.24 

288.18 

351.57 

291.39 


0.0016 

0.0019 

0.0023 

0.0029 

0.0039 

0.0058 

0.0059 

0.0103 

0.0112 

0.0116 

0.0101 

0.0068 

0.0038 

0.0020 

0.0012 

0.0009 


0.0009 

0.0010 

0.0010 


0.0008 

0.0014 

0.0017 

0.0027 

0.0034 

0.0030 

0.0016 


-0.49x10 

-0.27x10“ 

-0.27x10“ 


-3.99 

-3.66 

-3.33 

-3.00 

-2.66 

-2.33 

-2.00 

-1.67 

-1.34 

-1.01 

-0.67 

-0.34 

-0.01 

-0.32 

-0.65 

0.98 

1.32 

1.65 

1.98 

2.31 

2.64 

2.97 

3.31 

3.64 

3.97 
4.30 
4.63 
4.96 


i 


node 

time 

f 

u 

TABUS  VI 

(Fig.  12) 

fr  Q 

E 

010 

0 

564.94 

562.04 

563.49 

194.44 

0 

-2 

020 

3 

1137.75 

1133.55 

1135.65 

270.78 

0.77x10 

-2 

030 

6 

1699.31 

1694.34 

1696.82 

341.28 

0.28x10 

100 

9 

1137.85 

1133.56 

1135.71 

264.96 

0.48X10"4 

V^-lldB 

time 

f  V2 

V3 

V4 

P2/Pl  P3/Pl 

P4/Pl  FP 

12 

557.5 

-53.7 

-73.9 

13 

558.0 

-53.3 

-74.2 

14 

558.5 

-52.6 

-74.6 

15 

559.0 

-52.0 

-75.1 

16 

559.5 

-51.2 

-85.3 

17 

560.0 

-50.4 

-75.3 

18 

560.5 

-49.6 

-74.9 

19 

561.0 

-49.0 

-73.3 

20 

561.5 

-48.2 

-71.6 

21 

562.0 

-47.7 

-69.7 

22 

562.5 

-47.1 

-67.2 

23 

563.0 

-46.3 

-64.7 

24 

563.5 

-45.5 

-61.3 

25 

564.0 

-44.4 

-58.1 

26 

564.5 

-42.7 

-55.0 

27 

565.0 

-40.5 

-52.7 

28 

565.5 

-38.8 

-51.8 

29 

566.0 

-37.0 

-53.3 

30 

566.5 

-35.5 

-55.5 

31 

567.0 

-34.1 

-56.6 

32 

567.5 

-32.8 

-56.9 

33 

568.0 

-32.3 

-57.9 

34 

568.5 

-33.3 

-61.3 

35 

569.0 

-35.4 

-65.5 

36 

569.5 

-37.8 

-67.0 

37 

570.0 

-40.2 

-66.6 

38 

570.5 

-42.0 

-66.2 

39 

571.0 

-43.6 

-65.6 

40 

571.5 

-46.0 

-68.6 

-76.2 

-76.3 

0.0073 

0.0077 

0.0007 

-4.30 

-3.96 

-77.3 

0.0083 

-3.62 

-80.2 

0.0089 

-3.27 

-85.0 

0.0098 

-2.93 

0.0107 

C.0117 

-2.58 

-2.24 

-1.89 

-1.55 

-1.21 

0.0008 
a  nnna 

U«UX4Q 

0.0139 
n  A1AA 

0.0012 

UtUXHO 
n  m  co 

0.0015 
n  aaoi 

-0.86 

-0.52 

-0.17 

0.17 

0.52 

0.86 

U • UXDO 

0.0172 
0.0188 
0.0215 
n  aoaa 

a  nnn 

Um UUJi 

0.0044 

A  AAAI 

U. U^DU 

0.0335 

U  •  vvDJ 
0.0082 
A  AAQ9 

C 

"wO  •  J 

U.U4U/ 

U • UW74 

X  • 

-67.5 

0.0504 

0.0077 

0.0015 

1.55 

-67.0 

0.0596 

0.0060 

0.0016 

1.89 

-67.5 

0.0700 

0.0052 

0.0015 

2.24 

-63.3 

0.0813 

0.0051 

0.0024 

2.58 

-60.3 

0.0861 

0.0045 

0.0034 

2.93 

-59.3 

0.0072 

0.0031 

0.0038 

3.27 

-60.0 

0.0603 

0.0019 

0.0035 

3.61 

-59.7 

0.0457 

0.0016 

0.0037 

3.96 

-60.0 

0.0349 

0.0017 

0.0035 

4.30 

-60.3 

0.0282 

0.0017 

0.0034 

4.65 

-60.5 

0.0236 

0.0019 

0.0033 

4.99 

-62.0 

0.0178 

0.0013 

0.0028 

5.34 

mode 

time 

fu 

*L 

fr 

Q 

E 

010 

43 

565.47 

562.56 

564.01 

193.75 

0 

020 

46 

1138.89 

1134.67 

1136.78 

269.44 

0.78x10“' 

030 

49 

1700.94 

1695.60 

1698.27 

318.33 

0.37x10“' 

100 

52 

1139.00 

1134.76 

1136.88 

268.20 

0.81x10“ 
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TABLE  VII  (Fig,  13) 


mode 

time 

fu 

fL 

fr 

Q 

E 

010 

0 

561.56 

558.62 

560.09 

190.57 

0 

020 

3 

1137.45 

1133.15 

1135.30 

264.33 

0.14xl0_1 

030 

6 

1694.09 

1689.09 

1691.59 

338.39 

0.67xl0“2 

100 

9 

1137.58 

1133.26 

1135.42 

263.01 

O.llxlO-3 

Vj*-lldB  time 

f 

V2  V3  V4 

p2/Pi 

P3/Pl  P4/P; 

12 

554.1 

-53.6 

-74.7 

0.0074 

0.0007  - 

—  -4.21 

13 

554.6 

-53.0 

-74.7 

—  '  ■ 

0.0079 

0.0007  - 

—  -3.87 

14 

555.1 

-52.7 

-74.7 

0.0082 

0.0007  - 

—  -3.53 

15 

555.6 

-52.4 

-75.5 

0.0085 

0.0006  - 

—  -3.19 

16 

556.1 

-52.0 

-76.1 

. — 

0.0089 

0.0006  - 

—  -2.85 

17 

556.6 

-51.6 

-76.7 

0.0093 

0.0005  - 

—  -2.51 

18 

557.1 

-51.3 

-77.5 

—  .i 

0.0097 

0.0005  - 

—  -2.17 

19 

557.6 

-50.9 

-77.5 

0.0101 

0.0005  - 

—  -1.83 

20 

558.1 

-50.3 

-76.9 

—■ i  - 

0.0108 

0.0005  - 

—  -1.48 

21 

558.6 

-49.9 

-75.1 

0.0114 

0.0006  - 

—  -1.14 

22 

559.1 

-49.2 

-72.8 

0.0124 

0.0008  - 

—  -0.80 

23 

559.6 

-48.7 

-70.3 

0.0130 

0.0011  - 

—  -0.46 

24 

560.1 

-48.2 

-67.5 

0.0139 

0.0015  - 

—  -0.12 

25 

560.6 

-47.6 

-65.0 

— 

0.0148 

0.0020  - 

—  0.22 

26 

561.1 

-47.0 

-62.7 

0.0158 

0.0026  - 

—  0.56 

27 

561.6 

-46.3 

-60.4 

0.0172 

0.0034  - 

—  0.90 

28 

561.1 

-45.4 

-57.2 

0.191 

0.0049  - 

—  1.24 

29 

562.6 

-44.5 

-55.0 

. 

0.0211 

0.0063  - 

—  1.58 

30 

563.1 

-43.8 

-54.4 

0.0230 

0.0068  - 

—  1.92 

31 

563.6 

-43.1 

-55.0 

. . . 

0.0248 

0.0063  - 

—  2.26 

32 

564.1 

-42.4 

-57.5 

0.0269 

0.0047  - 

—  2.61 

33 

564.6 

-41.2 

-60.9 

0.0309 

0.0032  - 

—  2.95 

34 

565.1 

-40.0 

-64.5 

— 

0.0355 

0.0021  - 

—  3.29 

35 

565.6 

-38.6 

-67.9 

0.0417 

0.0014  - 

—  3.63 

36 

566.1 

-37.0 

-70.6 

0.0501 

0.0010  - 

—  3.97 

37 

566.6 

-35.2 

-73.7 

— 

0.0617 

0.0070  - 

—  4.31 

38 

567.1 

-33.2 

-71.0 

— ■ — 

0.0776 

0.0010  - 

—  4.65 

39 

567.6 

-31.5 

-66.3 

1 

0.0944 

0.0017  - 

—  4.99 

time 

fu 

*L 

fr 

Q 

E 

010 

43 

561.94 

559.01 

560.47 

191.35 

0 

020 

46 

1138.34 

1134.01 

1136.18 

261.97 

0.14x10' 

030 

49 

1695.35 

1690.34 

1692.85 

338.37 

0.68x10' 

100 

52 

1138.44 

1134.07 

1136.25 

260.07 

0.69x10' 
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TABLE  VIII  (Fig.  14) 


mode 

time 

f 

u 

fL 

fr 

Q 

E 

010 

0 

564.94 

561.90 

563.42 

185.40 

0 

020 

3 

1137.18 

1126.60 

1128.89 

246.64 

0.18xl0“2 

030 

6 

1703.17 

1698.34 

1700.75 

352.63 

0.62xl0“2 

100 

9 

1139.73 

1135.13 

1137.43 

247.16 

0.76xl0"2 

Vj—lldB 

time  f 

V2 

V3 

v4 

V*1 

p4/Pi  : 

12 

557.5 

-50.4 

13 

558.0 

-50.0 

14 

558.5 

-49.5 

15 

559.0 

-48.9 

16 

559.5 

-48.4 

17 

560.0 

-47.7 

18 

560.5 

-47.0 

19 

561.0 

-46.3 

20 

561.5 

-45.5 

21 

562.0 

-44.6 

22 

562.5 

-43.6 

23 

563.0 

-42.5 

24 

563.5 

-41.5 

25 

564.0 

-40.6 

26 

564.5 

-40.0 

27 

565.0 

-40.8 

28 

565.5 

-42.0 

29 

566.0 

-42.6 

30 

566.5 

-41.2 

31 

567.0 

-39.3 

32 

567.5 

-37.0 

33 

568.0 

-35.0 

34 

568.5 

-33.0 

35 

569.0 

-32.6 

36 

569.5 

-33.6 

37 

570.0 

-35.5 

38 

570.5 

-37.4 

-74.2 

— — - 

0.0107 

-74.2 

— —— 

0.0112 

-74.2 

— — 

0.0120 

-74.3 

— 

0.0127 

-74.5 

— 

0.0136 

-74.3 

0.0147 

-73.9 

0.0158 

-73.0 

0.0172 

-71.1 

0.0188 

-68.5 

— - 

0.0209 

-65.3 

0.0236 

-62.1 

— 

0.0267 

-59.3 

— 

0.0299 

-57.0 

-76.0 

0.0331 

-54.7 

-71.8 

0.0355 

-54.0 

-70.5 

0.0324 

-54.1 

-71.3 

0.0282 

-54.7 

-70.2 

0.0265 

-54.7 

-64.8 

0.0310 

-55.3 

-63.0 

0.0387 

-57.4 

-67.2 

0.05C1 

-58.7 

-72.3 

0.0631 

-58.0 

-75.5 

Q.0794 

-57.5 

-79.0 

0.0832 

-58.5 

0.0741 

059.9 

.  ■  . 

0.0596 

-63.5 

— — — 

0.0480 

-3.92 

-3.59 

-3.26 

-2.93 

-2.59 

-2.26 

-1.93 


0.0010 

-1.27 

0.0013 

-0.94 

0.0019 

-0.61 

0.0028 

-0.28 

0.0039 

0.05 

0.0050 

0.38 

0.0065 

0.0009 

0.71 

0.0071 

0.0011 

1.04 

0.0070 

0.0010 

1.37 

0.0066 

0.0011 

1.70 

0.0065 

0.0020 

2.03 

0.0061 

0.0025 

2.36 

0.0048 

0.0010 

2.69 

0.0041 

0.0009 

3.02 

0.0045 

3.35 

0.0047 

3.68 

0.0042 

4.02 

0.0036 

4.35 

0.0024 

4.68 

mode 

time 

fu 

fr 

Q 

E 

010 

0 

564.94 

561.93 

563.43 

187.00 

0 

020 

3 

1131.26 

1126.70 

1128 . 98 

247.64 

0.19x10' 

030 

6 

1703.35 

1698.48 

1700.92 

349.34 

0.63x10 

100 

9 

1139.88 

1135.27 

1137.58 

247.14 

0.76x10' 

010 

0 

561.20 

558.27 

559.73 

190.78  0  1 

020 

3 

1125.25 

1120.66 

1122.96 

244.76  0.31x10  *  \\ 

030 

6 

1700.26 

1695.60 

1697.93 

364.21  0.11x10  1  | 

100 

9 

1143.88 

1139.80 

1141.84 

280.07  0.17X10”1 

lldB 

time 

f 

V2 

V3  V4 

p2/Pi 

P3/Pl  P4/Pl  **  1 

553.7 

554.2 

554.7 

555.2 

555.7 

556.2 

556.7 

557.2 

558.2 

558.7 

559.2 

559.7 

560.2 

560.7 

561.2 

561.7 

562.2 

562.7 

563.2 

563.7 

564.2 

564.7 

565.2 

565.7 

566.2 

566.7 

567.2 


-52.7 

-52.4 

-51.8 

-51.4 

-51.0 

-50.5 

-49.7 

-49.1 

-47.8 

-47.0 

-46.0 

-45.0 

-43.6 

-42.6 

-41.0 

-40.5 

-41.0 

-43.2 

-46.0 

-49.0 

-51.7 

-53.6 

-54.4 

-53.7 

-52.4 

-50.6 

-49.7 


561.45 
1125.52 
1700.77 
1144 . 24 


-77.8 

-78.2 

-78.2 

-78.4 

-79.0 

-79.5 

-79.5 

-78.9 

-75.7 

-72.7 

-70.0 

-67.3 

-65.3 

-64.0 

-61.7 

-69.4 

-68.5 

-58.3 

-58.8 

-58.9 

-57.9 

-55.9 

-53.9 

-51.5 

-50.0 

-50.7 

-67.8 

*L 


558.48 

1120.92 

1696.08 

1140.17 


0.0082 

0.0086 

0.0091 

0.0095 

0.0100 

0.0107 

0.0116 

0.0124 

0.0145 

0.0158 


0.0005 


0.0178 
0.0200 
n  mv 

0.0011 
0.0015 
n  mi  q 

n  nocA 

n  nno? 

u  •  U  AD** 

n  iTJT  £ 

n  nooo 

u. uoxo 

n  nm  o 

w • UJjJ 

n  rm  £ 

n  nm  ** 

v • wJxw 

n  MAn 

n  nnA^ 

u • / 
n  m  tq 

n  ohai 

U • Ui / o 
n  ni  o£ 

n  nnAi 

U.  UJLAO 
n  nnao 

n  nnA£ 

-69.0 

U • UU74 

0.0074 

0.0057 

0.0013 

-64.5 

0.0068 

0.0072 

0.0021 

-62.4 

0.0073 

0.0094 

0.0027 

-61.9 

0.-085 

0.0085 

0.0029 

-60.1 

0.0105 

0.0105 

0.0035 

-56.0 

0.0116 

0.0115 

0.0056 

f 

Q 

I 

r 

-4.17 

-3.83 

-3.49 

-3.15 

-2.81 

-2.47 

-2.13 

-1.79 

-1.12 

-0.78 

-0.44 

-0.10 

0.24 

0.58 

0.92 

1.25 

1.59 

1.93 

2.27 

2.61 

2.95 

3.29 

3.63 

3.97 

4.30 

4.64 

4.98 


559.96 

1123.22 

1698.43 

1142.21 


108.60 

243.91 
362.29 

280.92 


0.29x10 
0.11x1 o' 
0.17x10' 
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